Microscopic Hamiltonian for Zn or Ni substituted high temperature cuprate 

superconductors 
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We have derived the effective low energy Hamiltonian for Zn or Ni substituted high-T c cuprates 
from microscopic three-band models consisting of the most relevant Cu or impurity 3d and O 2p 
orbitals. We find that both scattering potential and hopping integral induced by impurities have 
a finite range but decay very fast with distance from the impurity. The Zn scattering potential is 
very strong and attractive for electrons. The Ni scattering potential is much weaker than the Zn 
case, resulting from the hybridization between Ni ions and O holes. This profound difference is 
due to neither the electric charge nor d-level location, but rather because of the interplay between 
the valence state of the impurity and the strong correlation background, ft gives a natural account 
for the unusual effect of Ni and Zn on the reduction of superconducting transition temperature. 
The interlayer hopping of electrons is highly anisotropic and nonlocal, determined by the in-plane 
electronic structure. This leads to a quantum interference of states from different sites and affects 
strongly the scanning tunneling spectrum perpendicular to CuCh planes. 



Substitution of divalent transition metal Zn or Ni for 
Cu in Cu-0 plane has offered a fruitful and challeng- 
ing direction in exploring the nature of correlated elec- 
trons in high-T c cuprates J|-|,||-||. The effects of 
these impurities on superconducting properties are un- 
usual and in many ways opposite to those observed in 
conventional superconductors. In particular, the T c re- 
duction induced by non-magnetic Zn ions is almost three 
times larger than that by magnetic Ni ions ||. Also, the 
temperature dependence of the penetration depth shows 
a much stronger pair-breaking effect of Zn compared with 
Ni The actual situation is more complicated than it 
appears: the NMR and susceptibility measurements do 
show the existence of an induced local moment around 
Zn impurities 0. There exists a vast literature devoted 
to explore the impurity effects on d-wave superconduc- 
tors p-[T^] , particularly to interpret the scanning tunnel- 
ing microscopy (STM) experiments |^Jl^-p^]. However, 
the key issue, namely why Zn impurity is a very strong 
scatterer (close to the unitary limit), while Ni is a rel- 
atively weaker scatterer, is not well understood. The 
earlier efforts to address this issue using exact diagonal- 
ization of the t — J model on finite clusters [|l7| and the 
slave-particle treatments of the resonance valence bond 
(RVB) state H,E9| show the importance of the strong 
correlation background, but are unable to elucidate the 
physics involved to the full extent. 

In this paper, we present a theoretical derivation for 
the effective low-energy Hamiltonians for both Zn and Ni 
impurity systems. The aim is to find an accurate descrip- 
tion for the impurity scattering potential and to clarify 
a number of important issues associated with the scat- 
tering phase shifts. We also present a detailed analysis 
for the anisotropic interlayer hopping process and discuss 
its impact on the scanning tunneling spectra of high-T c 



superconductors. We start from the three-band Hubbard 
model for the Cu-0 plane, and study systematically the 
influence of the impurity on the Zhang-Rice singlet |2C| ] 
formation. In this work, only a single Zn or Ni impu- 
rity system is studied. However, it is straightforward to 
extend the derivation to a system of many Zn or Ni im- 
purities, provided the direct interaction among them is 
negligible. 

Let us first consider the case of Zn impurity. Zn 2+ has 
a 3d 10 configuration and no spin. Since the valence fluc- 
tuation of Zn 2+ is small, we can model Zn as an inert non- 
magnetic impurity with no coupling with surrounding O 
or Cu states. In high-T c oxides, the low energy physics is 
governed by one Cu 3d x 2_ y 2 and two O (2p x , 2p y ) bands 
and the doped holes are predominantly on O (2p Xl 2p y ) 
states. Thus to study the physical effect of Zn, we should 
start from the following three-band Hamiltonian 

H = H +H U (1) 

H = Ep^pjpi + ^ \£dd\di + Udd^diidl^dii 

l i^i 



Hi = - t pd (p\di 4- d\pA , 



where d\ = (4 v 4l 



creates a Cu 3d T 2_„2 hole, and 



Pi 



(PipPil 



creates an O 2p x hole if I — i ± x/2 or 



an O 2p v hole if I = i ± y/2. (See Fig. 1 of Ref. |2| 
for the definition of relative phases of id x 2_ y 2 and O 2p 
orbitals.) io is the impurity position. Without loss of 
generality, we set Ed = and e p > 0. The real system 
is in the limit [/,; > e p » t p d- An on-site Coulomb re- 
pulsion term for O holes is neglected since we are only 
interested in low doping cases and the probability for two 
holes to occupy one O site is very small. 
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The O 2p holes can be classified according to their hy- 
bridization with Cu 3d x 2_ y 2 states as bonding and non- 
bonding Wannier orbitals. These orbitals are defined 
on the Cu lattice. The bonding orbital is defined by 
a » = 7^ Y,k e lkRi Pk 1 [cos(k x /2)p x , k + cos(ky/2)p y>k ], 
where (3k — yJcos 2 (k x /2) + cos 2 (fc y /2) and N is the num- 
ber of Cu sites. p x ,k and p V;k are the Fourier trans- 
formation of P i+ ^j2 an d Pi + ^/2' respectively. The non- 
bonding orbital is orthogonal to a,i and defined by bi — 
E fc e lkR '(3^ 1 [cos(k y /2) Px , k - cos{k x /2)p ytk \. With 
these Wannier orbitals, the Hamiltonian (|l|) can be ex- 
pressed as 

i i=£i a 

H\ = -tpd ^2 u (i - j) (ojdi + 4 a i) i 

where u(i) — J^k Pke lkRi ■ The non-bonding orbitals 
do not couple to the other states and therefore may be 
neglected. 

In the limit Ud S> e p 3> t p d, we can treat H\ as a 
perturbation and use degenerate perturbation theory to 
project the Hamiltonian onto the Hilbert space spanned 
by the low-energy states of H a . Up to the second order 
of perturbation, we find that an effective Hamiltonian for 
the system is given by 

H 'zn = -^Ujaja-j +t 12 4x d «x'aL'««^ ( 2 ) 

where t 12 = t 1 + t 2 , h = t 2 pd /e p , t 2 = t 2 pd /(U d - e p ), 
tij = t 2 [4<5,;.j + — u(i — i)u(i — j)] (the second 
Kronecker symbol means i,j must be the nearest neigh- 
bours) and Hi = Y]j u(i — j) aj. The coupling between di 
and cij at different sites is much smaller than the on-site 
one. Thus, we can take the approximation <Zj ~ u (0) a^. 
This is also an approximation implicitly assumed in the 
work of Zhang and Rice j?o) . 

In H' , the coupling constant between O 2p and Cu 
id x 2_ y 2 states is of the order (t 2 + ti) u 2 (0). This is a 
large energy scale compared with the hopping integral for 
the bonding O holes. An important consequence of this 
strong coupling between O 2p and Cu 3d x 2_ y 2 holes, as 
first pointed out by Zhang and Rice (|o) , is that an O hole 
will form a local spin singlet with a Cu spin. The binding 
energy of the Zhang- Rice singlet is (t 2 + ti) u 2 (0). The 
energy difference between the Zhang- Rice singlet and the 
corresponding triplet is 2 (t 2 + t\) it 2 (0). Since this differ- 
ence is much larger than the bandwidth of the O holes, 
the triplet hole band can be neglected in the study of 
low energy excitations. Thus we can further project the 
Hamiltonian onto the subspace spanned by the Zhang- 
Rice singlets with unpaired Cu spins. The resulting ef- 
fective Hamiltonian is then given by 

H'zn = £ V ^)4d t - ]T tffd]di (3) 
i i^j 



subjected to the constraint d\di < 1 when i ^ {q. At 
the impurity site, di a = (di a ^di i) is defined from the 
creation operator of the bonding O hole: di ocr = —aa\ Q —. 
We should emphasize that di a is not an annihilation op- 
erator for the Zn 3d x 2_ y 2 state. At the Zn site, there is 
no constraint imposed on di . Thus two electrons with 
opposite spins can occupy this site simultaneously. This 
is consistent with the fact that there are two electrons in 
the 3d x 2_ y 2 state of the Zn 2+ ion, although our starting 
Hamiltonian ( [j]) does not contain that state at all. More- 
over, if the total number of doped holes in the system is 
Nh, then it can be shown that d]d% — (N — Nh) + I. 
Thus the total number of conducting electrons is more 
than the nominal charge of the system by one. This im- 
plies that Zn 2+ can be effectively taken as an ion with 
one more charge than Cu 2+ , although both are divalent. 
It further suggests that we only need to count the dif- 
ference between the number of 3d electrons in Zn 2+ and 
Cu 2+ when the Friedel sum rule of phase shifts is used. 
In (||), the hopping integral is defined by 

,Zn _ tj±s: , «>oj r , hpi r (a\ 

L ij — 2 ° 4 Aoj'Ao ' \/2 

where tij — t 2 S^j^ — t 2 u(io — i)u(io — j)- The first term 
in tij is the usual hopping integral in a homogeneous sys- 
tem, and the second term is a small correction to the 
hopping integral by the Zn impurity. This impurity in- 
duced hopping integral is non-local but decays very fast 
with distance from the impurity. 

The impurity potential Vz n {i) is given by 
V Zn (i) = -hu 2 (i - i)5^ io - (2< 2 + ti) u 2 (0) 5 itio . (5) 

This is an attractive potential. It is certainly not of a 
(^-function form, but away from the Zn impurity, Vz n (i) 
decays very fast, roughly as \i — io|~ 6 . At the impurity 
site, |Vz„(io)| ~ 3.67(ii + 2£ 2 ) is about two orders of 
magnitude larger than the potential at the four nearest 
neighboring sites of Zn, |Vz n (i)| ~ 0.0785£2- It is also 
more than one order of magnitude larger than the hop- 
ping constant £2/2, in consistent with the STM measure- 
ment data PJl5[| . This shows clearly that Zn impurity is 
indeed a strong potential scatterer. 

The strong attractive potential induced by Zn is in 
fact due to a strong repulsion of Zn 2+ to the bonding 
O holes at the impurity site. This strong repulsion of 
Zn 2+ to O holes results from two distinguished physical 
effects. First, since Zn 2+ has no spin, O holes cannot gain 
energy by forming a Zhang-Rice singlet with a 3d x 2_ y 2 
spin at the Zn site, this costs the O hole an energy of 
(t 2 + ti)u 2 (0). Second, since there is no hopping between 
O 2p and Zn 3d states, the O hole also suffers from a ki- 
netic energy loss of t 2 u 2 (0). The total energy loss of an 
O hole at the impurity site is thus (2t 2 + ti)u 2 (0). This 
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serves effectively as a repulsive potential to the O holes, 
or equivalently an attractive potential to the di elec- 
trons, at the Zn site. This picture is completely different 
from the naive interpretation for the origin of the attrac- 
tive potential of Zn in literature, which claims that the 
potential is attractive because Zn 2+ has one more elec- 
tron in its 3d x 2_ y 2 orbital than Cu 2+ and the 3d energy 
levels of Zn 2+ are lower than those of Cu 2+ . We would 
like to particularly point out an intriguing connection be- 
tween the scattering potential of Zn and the physics un- 
derlying the Zhang- Rice singlet. It suggests that through 
the measurement of the Zn scattering potential, we can 
even gain direct information on the binding energy of the 
Zhang- Rice singlet. 

The Hamiltonian (^) contains the most important in- 
formation about the Zn scattering potential. However, to 
obtain the full effective Hamiltonian, we should include 
the fourth order of perturbation as in the derivation of 
the standard t — J model. In particular, two terms which 
are absent in the second order of perturbation should 
be included. One is the exchange energy between two 
Cu 2+ spins on nearest neighboring sites. The other is 
the hopping of electrons between next or next-next near- 
est neighboring sites. This term has strong effect on the 
structure of Fermi surface. There is also a three-site hop- 
ping term, which is generally ignored for simplicity. The 
other terms generated from the fourth order of pertur- 
bation are small corrections to the terms in (|^) and can 
be omitted. Thus the full effective Hamiltonian for the 
Zn-substituted system is 

H Zn = H% n -j2 tu4dj + E JS *- s i> ( 6 ) 

where Si is the Cu spin operator, and t'^ = t' if 
are next neighbors, or t" if (i,j) are next-next neighbors, 
or zero otherwise, t' and t" always have opposite signs. 
They are of the same order as the exchange energy J. 

Now let us turn to the Ni case. As for Zn 2+ , Ni 2+ is 
also in a divalent state. It has 8 electrons in its 3d states. 
In particular, it has one electron in the 3d 3z 2_ r 2 state 
and one in the 3d x 2_ y 2 state. These two electrons form 
a spin triplet via a strong Hund's rule coupling. Like 
Zn 2+ , Ni 2+ is also very stable. However, unlike Zn + , 
the hybridization between the Ni (3d x 2_ y 2, 3d 3z 2_ r 2j and 
O (2p x , 2p y ) states is strong since there is only one elec- 
tron in each of these 3d states. 

The appropriate three-band model for the Ni system 
is defined by the Hamiltonian 

H m = H + Y, c+Ca ~ E C (pW + h.c.) 

a (U )a 

-JHc\-d • 4|c 2 + E U a lc U C <xt c ll C <xl: ( 7 ) 
a 

where H is the same as defined in ([j]) . a = 1 or 2 denotes 
the Ni 2+ 3d x 2_ y 2 or 3d 3z 2_ r 2 state. c a = (c a f, c a |) is the 
electron operator for these Ni 2+ states. e„ l , U^ 1 , and 



t^ 1 are the energy level, the on-site Coulomb repulsion, 
and the hopping integral of the corresponding 3c? state 
of Ni 2+ , respectively. Jh is the Hund's rule coupling 
constant. 

In (|t]), t^ 1 is much smaller than Jh- Thus this t% % 
term can be also treated as a perturbation. The effective 
one-band Hamiltonian for the Ni impurity can be derived 
following the same procedure as for the Zn impurity. The 
difference is that an O hole can now form a local spin 
doublet (an analogue to the Zhang-Rice singlet) with the 
Ni 2+ spin at the impurity site. This difference is quite 
important, and it leads to a relatively weaker scattering 
potential for Ni. The effective one-band Hamiltonian we 
obtained for the Ni impurity is given by 

+ E Vm(i)d\di + J ij s i ' s i ( 8 ) 

where T = (Ti,To,T_i) is the annihilation operator 
for the three spin triplet states of Ni 2+ , and di = 
(di i , di i) is the annihilation operator of the local spin 
doublet formed by the Ni 2+ spin and the bonding O 
hole at the impurity site. As for the Zn case, di 
is not an electron operator for the 3d state of Ni 2+ . 

T} = (4 a Av(<Al+<l<)/^<Ai)- When 
i 7^ io, the constraint d\di < 1 is imposed. At the im- 
purity site, the constraint becomes d\ o di + T^T = 1, 
namely the impurity site must be either in a spin triplet 
state of pure Ni 2+ or in a Zhang- Rice-type spin doublet 
state. The second term in H^i describes a swap process 
between the spin triplet and the spin doublet at the Ni 
site by absorbing or releasing an electron at the other 
site. The exchange constant Jij = J if both i and j are 
not equal to io- For the bonds connecting the Ni site 
is different from but of the same energy scale as J. 
The hopping integral tf^ is 

tf/ = t% + {U - t 2 ) u{i Q - »)«(t - j), (9) 

where is the hopping integral as in the standard t—J 
model without impurities, including the hopping between 
nearest, next nearest and next-next nearest neighbors. 
The second term in tf^ 1 is the hopping integral induced 
by the Ni impurity. This term is smaller than the corre- 
sponding term in the Zn case. 

The scattering potential of Ni reads 

V m {i) = (2t 2 +ti- ~t 3 - ~* 4 ) u 2 (0)S l:to 

- (*2 - t&) u 2 (i Q - i)Si^ io , (10) 

where t 3 = £ Q (t£*) 2 /(e p - £ f + J H /i) and t 4 = 

J2 a (*a f ') 2 /(^ ~ £ p + e - 1 + J W4). The t 3 term origi- 
nates from the hopping between Ni 2+ 3d and O 2p states. 
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The ti term comes from the energy gain by the forma- 
tion of the local spin doublet at the impurity site. These 
terms compensate or even overcome the energy loss re- 
sulted from the t\ and t% terms. Thus the scattering 
potential of Ni is much weaker than Zn. 

The total number of unpaired electrons is now deter- 
mined by the formula £ i#io d\di = N - N hole - T^T. 
If the scattering potential is strong and repulsive to the 
spin doublet at the Ni site, then this site is essentially 
in the spin triplet state, i.e. T^T = 1, and the conduct- 
ing electron number is one less than the total number 
of charge. In this case the Friedel sum rule should be 
used as if one charge is absorbed by the Ni 2+ ion. If, 
on the other hand, scattering potential is strong but at- 
tractive to the spin doublet at the Ni site, then this site 
is predominantly in the spin doublet state and the total 
number of conducting electrons is equal to the total num- 
ber of charge. In this case, the scattering phase shifts of 
high angular momentum channels become important and 
the commonly adopted T-matrix approximation 22 is no 
longer applicable. In real systems, the Ni site might be 
in a mixed valence state, with finite probabilities in both 
spin triplet and doublet states. In this CcLSG, el thorough 
analysis of the Ni scattering becomes very complicated, 
even without considering its magnetic coupling with Cu 
spins. 

By comparing the results for Zn and Ni cases we clearly 
see the difference: Zn is a very strong scatterer, while Ni 
is a weak one. This profound difference is due to nei- 
ther the electric charge nor <i-level location, but rather 
because of the interplay between the valence state of the 
impurity and the strong correlation background. The 
underlying physics is best revealed using the three band 
model and the concept of Zhang- Rice singlet. The pene- 
tration depth measurements give us strong hints that Zn 
materializes a unitary scattering, whereas Ni is a scat- 
terer close to the Born scattering limit ||. The present 
microscopic derivation provides a very good justification 
for this scenario. We did not calculate explicitly the mag- 
netic moment distribution and the contribution of spin- 
flip scattering around a Zn or Ni impurity. However, 
we believe that these problems can be addressed by self- 
consistently solving the effective Hamiltonians for Zn or 
Ni impurities . The fact that the suppression of T c scales 
universally with the increase of the residual resistance [^) , 
and the rough estimate of the dominance of the poten- 
tial scattering over magnetic scattering 0, do support 
that the spin-flip scattering induced by Zn or Ni impu- 
rities is a secondary effect compared with the potential 
scattering, especially for the Zn case. 

Finally, let us discuss the problem about the inter- 
layer hopping of electron around an impurity. This prob- 
lem is of great importance for the understanding of vari- 
ous experimental results related to impurities, especially 
the scanning tunneling data in the superconducting state 
around Zn or Ni ||. In high-T c materials, the interlayer 
hopping is fulfilled through a virtual hopping process as- 



sisted by Cu 4s orbitals. That is an O hole first hops 
to a Cu 4s orbital and then hops to other orbitals out 
of the CuC>2 plane (2^j2^|. At the impurity site, the in- 
terlayer hopping is assisted by the Zn 4s or Ni 4s state. 
Since the 4s orbital is orthogonal to the 3d x 2_ y 2 orbital 
at the same site, there is no hopping between these two 
orbitals. This 4s-orbital assisted hopping has profound 
consequences on c-axis dynamic properties [p3| , pH . 

The 4s-orbital assisted hopping in each Cu02 layer is 
described by the following Hamiltonian: 



s \pi 



pjs, 



(11) 



where (-)"" = 1 if I = i ± x/2 and {-) a " = -1 if 
I = i ± y/2. sj = (si., st.) creates a hole in the local 4s 
state. tf s is the hopping integral between a 4s and its 
neighboring O 2p orbitals. By projecting this Hamilto- 
nian onto the low energy subspace containing Zhang- Rice 
singlets, we find that H c can be expressed as 



H C = J2 ***** A + h - 



(12) 



where s\ = (— 3^, s,-f) is the electron creation operator 
for the 4s state at site i, 



D, 



l 

V2 



for the Zn impurity model, and 



V27\, To 

V2 UJ V3 Ul ° V To, 



(13) 



(14) 



for the Ni impurity model. In Eqs. ( |l3| ) and (|l4|), 

F id = ^-^2p^(cosk x - cos k y )e lk ^l (15) 



It is straightforward to show that Fij — when i = j 
or when i — j is along two diagonal directions. Thus in 
the effective one-band model, an electron in the CuC>2 
plane cannot hop to a 4s orbital at the same site. This 
is a peculiar but important feature of the c-axis hopping 
integral. 

The local tunneling current perpendicular to Cu02 
planes as measured by the STM is proportional to the 
probability of electrons in the Cu02 plane hopping to 
the local 4s orbital from the low energy conducting band. 
From Eq. (12) one can show that the tunneling conduc- 
tance at site i is determined by the following non-local 
matrix element 



9i (V) cx (Di5(eV 



(16) 



where V is the applied bias voltage and Hi mp is the sys- 
tem Hamiltonian as given in Eqs. M) or (0). This formula 
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shows clearly that the scanning tunneling conductance 
is not simply a measure of the local density of states in 
high-T c superconductors. In particular, as Fij = when 
i = j, there is no contribution to gi(V) from the low en- 
ergy excitations at site i. This is a peculiar but intrinsic 
property of the in-plane electronic structure. It is not 
due to the blocking effect as suggested in Ref. [^6|. In 
Eq. (g) or (|l^) , if j in the summation only takes the 
four nearest neighboring sites of i, we then obtain the 
result given by Martin et al [jL5[ . Their result is a good 
approximation to Eq. (|l^) since Fjj decays very fast 
with the distance \j — i\. Eq. ( )l6| ) is valid independent 
on the structure of the insulating layers between the top 
Cu0 2 plane and the STM tip. Thus the normalized STM 
spectrum should be the same no matter whether the tip 
is above an insulating layer or directly above a Cu02 
layer. 

In conclusion, we have derived the effective low en- 
ergy Hamiltonian with the impurity scattering potential 
for Zn or Ni substituted high-T c cuprates. Our results 
reveal an intrinsic connection of the impurity scattering 
potential with the correlation effect of electrons and set 
a solid microscopic starting point towards thorough un- 
derstanding of the impurity effects in these materials. A 
number of observable physical quantities can be calcu- 
lated explicitly based on this microscopic model. The 
STM spectrum perpendicular to CuC>2 is not simply a 
probe of the local density of states since the interlayer 
hopping is intrinsically anisotropic and non-local. 
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